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ABSTRACT 

We provide progenitor models for electron capture supernovae (ECSNe) with 
detailed evolutionary calculation. We include minor electron capture nuclei us- 
ing a large nuclear reaction network with updated reaction rates. For electron 
captures, the Coulomb correction on the rates is treated and contribution of 
neutron-rich nuclei are taken into account in the nuclear statistical equilibrium 
(NSE) composition. We calculate the evolution of the most massive super asymp- 
totic giant branch stars and show that these stars undergo off-center carbon 
burnings and form ONe cores at the center. These cores get heavier up to the 
critical mass of 1.367 and keep contracting even after the initiation of 0-|-Ne 
deflagration. Though inclusion of minor electron capture nuclei causes convective 
URCA process at the contraction phase, such process will have minor effect on 
the evolution. On the other hand, electron captures by neutron-rich nuclei in the 
NSE region have more significant effect. Also, we discuss the uniqueness of the 
critical core mass for ECSNe and the effect of mass loss on the plausibility of our 
models for ECSN progenitors. 

Subject headings: nuclear reactions, nucleosynthesis, abundances, — stars: evo- 
lution — stars: interiors — supernovae: general 



1. Introduction 

Electron capture supernova (ECSN) is a quite distinct class in core collapse supernova 
(CCSN). While an ordinary CCSN is considered to be triggered by a collapse of an Fe 
core, which is located at the center of the star with an onion-like compositional structure, 
a progenitor of an ECSN is a super asymptotic giant branch (SAGB) star that has a core 
mainly co nsists of oxyge n and neon, surrounded by a thin helium shell and diffuse hydrogen 



envelope (iNomotd 119871 ). In an ONe Chandrasekhar mass core, electron capture reactions 
on ^^Mg and ^°Ne heat the surroundings. As a result, O and Ne burning ignites at the center 
and generates energy, and 0-|-Ne deflagration propagates outward. However, the released 
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energy is too small to explode the highly bound core (IMiyaji et al.lllQSd ). Further electron 
capture reactions in the central NSE region (Nuclear Statistical Equilibrium) accelerate core 
contraction . Finally, a proto-n eutron star forms and the star ends up as a weak type II 
supernova ( iKitaura et al.l 120061 ) . 



The most distinct point of the progenitor may be its contrasting structure of a highly 
conce ntrated core and a diffuse envelop e. With a pro g enitor model of an ECSN (INomoto 
19871 ). a hydrodynamical simulation by Kit aura et al.l (120061 ) showed that the delayed ex- 
plosion powered by neutrino heating takes place even in a one-dimensional calculation, 
and the successful explosion is confirmed by recent multi-dimensional calculations as well 
( Janka et al. 2012[). From these self-cons istent calculations, properties of ECSNe such as 
nucleosynthesis ( IWanajo et al.l |20 111 l2013l ) have been studied. 



Not only for a theoretical aspect, but also for observatioi is, a model of ECSN has impor- 
tant suggestions. Sorn e low luminosity SNe, e.g., SN1997D ( iTuratto et al.l 119981 ) . SN2005cs 
( jPastorello et al.ll2006l . |2009| ). can be explained by the explosion model of an ECSN which 
has a low explosion energy and synthesizes small amount of ^^Ni. Also observed peculiar 
compositions in the well-known Crab nebula, such as abundant He and less abundant C 
and O, indicate tha t the Crab supernova SN1054 arose from a collapse of a SAGB star. 
( INomoto et al.lll982[ ). Type Iln SN, which is a SN explosion enshrouded by a dense circum- 
stellar medium, can be explained by an ECSN as well as a CCSN from very massive star 
that has experienced an intense mass loss phase. Re cently, a progenitor of a dust-enshrouded 
transient SN2008S is found in a pre-explosion image (iBotticella et al.ll2009l ) and it would have 
a mass of ~10 M©, which is a plausible mass for an ECSN progenitor. 

However, there has not been a consistent progenitor calculation from zero age main se- 
quence (ZAMS) to collapse because of the numerical difficulties to calculate the full evolution 
of SAGB stars. The main difficulties are off-center C burning, thermal pulses, contraction of 
a highly degenerate core, calculation of electron captures, and propagation of deflagration. 
These phases had been separately studied by several authors. 

Th e theoretical work on a collapsing ONe core was initiated by Nona o to and collaborators 
in 80 's dMivaii et alll98ol : iMivaii fc Nomotolll987l : lNom"otolll984j . E987[ ). iMivaii et all Jl98o[ ) 
investigated the effect of electron capture on ^'^Mg and ^°Ne and showed that these effects 
can be summarized as follows: Firstly, reduction of the electron mole fraction induces core 
contraction. Secondly, reduction of the electron mole fraction reduces the Chandrasekhar 
mass. And finally, electro n capture affect s the energy equation by both endothermic and 
exothermic way. Moreover, iNomotd ( 119871 ) followed the core evolution either after the initi- 
ation of 0-|-Ne deflagration using a He star model and provided the progenitor model for an 
ECSN. This model has been the only one that can be used in an explosion simulation. In 90's, 
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non-explosive evolutionar y calculation of solar-metal SAGB stars was investigated by Garcia- 



Berro and collaborators ( Garcia-Berro &: IbenI Il994j : iRitossa et al 



1996 



Garcia-Berro et al. 



I997I : Ilbenlll997l : iRitossa et aljll999h . Non-solar metallicitv (|Gil-Pons et ali boosl: ISiesslboOTh 



as well a s detailed p hysics such as overshooting (iGil-Pons et al.ll2007l ) and thermohaline con- 
vection (jSiessI 120091 ) are considered in recent studies. An off-center carbon burnin g and 



thermal pul ses, which require a high cost calculation, were extensively investigated by 



Siess 



( I2OO6I . I2OIOI ). The contraction phase of an ONe c ore was investigated b y recent works as 
well, concerning t he different nuclear rea ction rates (IHashimoto et al.lll993r). different convec - 



tive assumptions (iGutierrez et al.lll996l ). and different compositions ( iGutierrez et al.ll2005l ). 
These simulations stopped at the ignition of O burning and the continuous deflagration phase 
had not been calculated. 

Thus the main purpose of this work is to calculate a progenitor model for an ECSN 
from a detailed stellar evolutionary si mulation . The calculation treats its main sequence 
phase, which is omitted in the work by iNomotd (119871 ). Also an off-center C burning phase, 
which is one of the important improvement in the evolutional theory for SAGB stars, is 
followed in detail. Using a large nuclear reaction network, we include minor nuclei that are 
synthesized during C burning phase and additional electron capture reactions on these nuclei 
are treated with the Coulomb correction on rates. Update d nuclear reaction r ates, especially 
the new electron capture rate for the NSE composition by lJuodagalvis et al.l ( I2OIOI ). is taken 
into account. The growth of an ONe core mass is assumed to be a result of a stationary 
helium burning. This enables us to avoid numerical difficulties on shell He burning and to 
investigate the full evolution of an ONe core. 

The critical core mass for ECSNe, with which an ONe core can no longer exists as 
a stable object because of the initiation of electro n captur e react ions, is considered to be 
a uniquely determined quantity an d the value by iNomotd (119871) is u s ed for estirnates o f 
the initial mass range for ECSNe (Siess 2007 : Poelarends et al.l 2008 : Pumo et al. 2009). 
However, the uniqueness on various parameter settings has not been confirmed yet owing to 
lack of calculations. Moreover, update of both numerical prescriptions and physical effects 
such as diffusive convective mixing and the Coulomb correction on electron capture rates 
possibly affect the value. Thus investigation on the parameter dependence of the critical 
core mass for ECSNe is intended in this work as well. 

We report stellar evolution of the most massive SAGB stars with solar composition from 
its main sequence through ONe core contraction, and further deflagration phase in which 
the central NSE region extends outward. In the next section, the method of calculation and 
input physics are explained. The evolutional results from ZAMS phase to the formation of 
the ONe core are presented in § 3.1., and core contraction including the propagation of NSE 
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front are explained in § 3.2. Discussions and conclusions are given in § 4. 



Methods 



The code was modified from lYoshida &: Umedal (1201 ll ) (hereafter, referred as Y&Ull) 
in order to calculate the late phase of ONe core evolution. The mixing length was set to 
be 1.5 times pressur e scale hight. Other assumptions and description of the code is given in 
Umeda et al.l (120 12[ ). In the foUowings, we explain other modified points. 



2.1. Spacial resolution 

Mesh points are automatically replaced to achieve required resolutions for calculations 
of stellar structures. Different conditions are set for different environments as well as different 
evolutional phases. Especially, for the case of propagation of shell C burning, which requires 
careful treatment of mesh refining, 

(i) |AlogP| < 0.05 

(ii) |AlogT| < 0.1 
(in) |Alogr| < 0.1 

(iv) V(AL,)VL2 < 0.1 

(v) AMr/M < 10"^ 

are taken as the constraints, where A/ means the difference of / between two mesh points 
and symbols have their usual meanings. 



2.2. Nuclear reaction netvi^ork 



We include 300 species of nuclei in a reaction network from n, p to Br. That is used 
for the calculation of distribution evolution of chemical composition and also nuclear energy 
generation. Newly added nuclei from Y&Ull (see also lYoshida et al.l (120131 )) are shown in 

Fig. m 

When the temperature exceeds the value of 5.0 x 10^ K, the composition in the region 
is assumed to achieve the NSE composition and the composition is used for the calculation 
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of the thermodynamical quantities and nuclear generation rate. In a NSE region in a highly 
dense ONe core, both density and electron mole fraction changes its value in wide ranges 
of (9.0 < log p < 11.0) and (0.25 ^Y^. < 0.5). In order to calculate the NSE composition 
consistently in such a wide parameter space, 3091 nuclei are treated for calculations of NSE 
ab undance (Fig. El). For ele ctron capture reactions in the NSE region, we applied the rate 
by lJuodagalvis et al.l (120101 ) which takes into account roughly 2700 nuclei including very 
neutron-rich and heavy ones with screening corrections. Though the contributions from 
positron capture and decay as well as /3~-decay in the NSE region may affect the stellar 
evolution, we omit these effects because of the lack of data tables. 



2.3. Convective criterion and diffusive mixing approximation 



In our code, convective mixing is treated as a diffusive process, and mixing beyond 
the convective boundaries such as overshooting is not treated. Convective boundar ies are 
determined by the Schwarzschild criterion, and semi-convective diffusion coefficient by lSpruit 
fll992h is also applied, thus 



fscD 



\ulcv for Wrad " ^ ad > niin(0, |V^) 

^--d-^-'i for < Vrad - V„ 



therm 



5 



ad — 5 * M' 



where u and lev are convective velocity and scale length determined by the mixing length 
theory, respectively. Dthp rm. = ^ '^'^^'^J^ is the thermal diffusivity, and fsc is a free parameter 



taken as 0.3 according to lUmeda fc Nomotd (120081 ). In order to take into account the effect 
of degeneracy of electrons both into the criterion and the coefficient of convective mixing, 
we define 







by extending the work in iKatd (119661 ) , where 

c^lnp _ dlnp _ dlnp 



5 



dlnT' 



dlnpi 



dlnpi 
d\nP 



dlnpe 
(MP 



(2) 



(3) 



Note that thermohaline convection which would take place in a region of f < is not 
treated for the sake of simplicity. 
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2.4. Electron capture and /3 -decay 



2.4- 1- The energy equation 

In a dense ONe core with density > 10^ g/cm^, the energy release by electron capture and 
/3~-decay of C burning products becomes important in the energy equation. The divergence 
of the energy flux Lr in the stellar equation is written as 

dLr 



dMr 



e„ - e,. 



-T 



E dsk \ 
. dt J 



(4) 
(5) 



where Sk is the s pecific entropy for th e k-th particle: nuclei, electrons, and photons ( iMiyaji et al. 



(119801), see also 



Ritossa et al.l ( Il999l )). e„ and e^, are the nuclear energy generation rate and 



the neutrino energy loss rate due to processes other than electron capture and f3 -decay 
respectively, e^eak is the energy generation rate of electron ca pture and /j~-decay, an d emix 



repre sents a cooling term owing to the work by convection ( ICouch fc ArnettI Il975l : llben 
19781 ). For radiation, the chemical potential is assumed to vanish, while the one of nuclei 
is ignored up to the achievement of NSE. For electrons, the chemical potential affects the 
energy equation through both e^eak and emix- Here we omit the description of other weak 
processes such as positron capture and /3"'"-decay because of small effects of these reactions. 

The energy generation rate of both electron capture and /3~-decay consists of three 
terms: mass difference, ne utrino emission, an d chemical potential excluding subatomic en- 
ergy of relevant particles ( IMiyaji et al.lll980[ ). Since we ignore the chemical potential of 
nuclei in the energy term, the total rate becomes 



^-weak 



dt 



(6) 



where Arrij is the mass difference between parent and daughter nuclei, (E^^j) is the mean 
energy of emitted neutrino, /ig is the chemical potential of electron, and 



e,weak.j 

~~dt 



is a time 



derivative of a specific electron number density owing to the j-th electron capture or 
decay, respectively. The sign of the neutrino energy is taken positive for electron capture 
and is taken negative for /3~-decay. The specific electron number density changes according 
to 



dt 



(7) 



where Uj denotes the specific number density of the j-th nucleus, and the sign of the reaction 
rate of electron capture or /3~-decay, A^, is taken negative for electron capture and is taken 
positive for /3~-decay, so that the emission of neutrino should always take some energy away. 
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According to llbeni (Il978l ) , we define 



Fe{Mr 



or 

where Fe{Mr) is the flux of electrons and is defined as 

Fe{Mr) = (47rrV)/^™x^. (9) 

Practically, the electron flux is calculated explicitly in our code, integrating the result of the 
chemical mixing, thus, 

(47rpr^)F, = ^^dM, (10) 
Jo at 

where is a time derivative of a specific electron number density owing to chemical 

mixing. In this expression, the energy loss vanishes at the boundaries of the convective 
region where the net flow of electrons should be zero. 



2.4-2. Correction for the reaction rate of electron capture and (3 -decay 



The Coulomb screening fo r the electron capture rate (ICouch fc Loumoslll974j : [Gutierrez et al 



19961 : lJuodagalvis et al.ll2010[ ) is taken into account in our code. Originally, the rate of the 



electron capture of j-th nucleus \ec,j is written as 



1 



E 



states 



1 



1 + exp 



. kT ' 



-dep 



where ae.r.,i is the electro n capture cross section and is the threshold energy of the reaction 
( lJuodagalvis et al.ll2010l ). This equation imphes that Aecj rapidly increases when /Xe exceeds 



e°, since electrons obey the Fermi-Dirac statistics. In highly dense regions, the equation 
of state of nuclei should include the Coulomb correction, and this affects the pressure, the 
entropy and also the chemical potential owing to the change of the Helmholtz energy. Because 
the chemical potential depends on a proton number Z, change in Z through the electron 
capture (A, Z) + e~ — )■ (A, Z — 1) + v^. will change the threshold energy of the reaction by 
the amount 



ec 

/i(Z) 



.z 



r 



0.9 + c,(- 



1/3 



'ZJ 



+ 



«o + "1 ( ^ 



(12) 
(13) 
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where Tz = Z^/^Z^/^^^^, = {j^Y^^ is the Coulomb couphng parameter for Z nucleus, Z 
is the mean charge of ions, a nd the four constants are ci = 0.2843, C2 = —0.054, do = —9/16, 



and di = 0.460 respectively (jPeWitt et al.l 119731 ). Therefore the effective threshold energy 
becomes 

4e// = d + AeL- (14) 

We assume that the energy distribution function for degenerate electrons does not change 
its shape in a small difference in electron density. Thus to take into account the effect of 
the Coulomb correction, we first evaluate the effective electron density that reproduces the 
effective chemical potential as the corrected Fermi energy, 

eF,eff = eF-Ae°^ (15) 

= f^e,eff{{pYe)eff)- (16) 

Then the reaction rate with effective electron density and fixed temperature is applied as 
the corrected electron capture rate, 

Aecj,e// = Xec,j{{pYe)eff,T). (17) 

Since Aeg^ is positive, this correction increases the effective threshold energy, reduces the 
effective Fermi energy, and thus reduces the rate of electron capture. 

The same correction is also apphed to /3~-decay. In the case of {A, Z) — )■ (A, Z + 1) + 
e~ + Vei change in the threshold Ae^^ is given as 

Ae°, = ^^{Z + l)-^,{Z) (18) 

which takes negative value and reduces the threshold energy. For /3~-decay, the decay rate is 
correlated to the number of unfilled states of electron in which the kinetic energy is between 
the threshold energy and the chemical potential, i.e., pe < < ^bd,eff- Therefore, the 
decrease of the effective threshold energy reduces the reaction rate of /3~-decay as well as 
electron capture. 



2.5. An approximate treatment of core growth 

After the completion of carbon burnings in a forming ONe core, we assumed a constant 
core mass growth in order to avoid high cost calculation of shell helium burning, and we 
assumed that the envelope will remain and shell He burning will continue until core collapse. 
The radiative-zero approximation was applied to the edge of the core as boundary conditions. 
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Under this assumption, three rates, 1.0x10"^ Mq/jt, 1.0x10"^ Mq/jt, and l.Ox 10"'' Mo/yr 
are taken as the core growth rate. The middle one of l.Ox 10~^ Mq/jt will be the most typical 
rate as the core growth from shell helium burning. This rate is consistent with the work by 
Nomoto (1987) in which the steady helium burning is assumed, and also with recent study 
by Siess (2010) and Poelarends et al. (2008) in which the TP phase is calculated. The results 
shown in § 3 are taken from the case of this typical rate. 



2.6. Late phase of core evolution 

When the time scale of co re evolution becomes shorter than that of convection, the well- 
known mixing length theory ( Bohm-Vitensel E-9581 ) which assumes a stationary convection 



becomes invalid. In order to determine the te mpera t ure gr adient in such cases, the time 



dependent mixing length theory formulated by lUnnd (119671 ) is adopted in our calculation. 
In this scheme, two differential equations on time for convective velocity Vcv and temperature 
fluctuation AT are given as 

/d_ ^ AT / dhgp \ dP_ 

Vcit^Z,,/2r™ 2pT\d\ogT)pdr' ^ ' 

^ + -^)aT = ^f-i£^)T(V- Vad) (20) 



where V — Vad repr esents the exce ss of temperature gradient compared with the adiabatic 



gradient. Following iNomotd (119841 ). we take the length scale of time dependent mixing l^v 
as the smaller one among the radial distance and the local mixing length. The convective 
energy flux Fcv is written as 

= cppATu (21) 

where cp is the specific heat at constant pressure. Then identical to the mixing length theory, 
equations of total luminosity 

Lr = Lrad + ^Trr^F^^, (22) 
Uad = ^;7^v (23) 

are solved to get the temperature gradient. 

As the timescale of evolution decreases and becomes comparable with the free-fall 
timescale, the assumption of hydrostatic structure becomes invalid. In this work, an inertia 
term is included in the equation of motion and also in the radiative temperature gradient 
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flHeger et al.lboooh as 



dP 

dMr 

^ rad 



3kPL 



1 + 



d r 



IGnacGMrT^ 



1 + 



-1 



GMr df^ 



(24) 
(25) 



3. Results 



We calculated the evolution of 10.4-11.2 Mq stars with a metallicity of Z=0.02, from 
ZAMS to 0+Ne deflagration for 10.4-10.8 Mq models and from ZAMS to off-center Ne 
ignition for 11.0, 11.2 Mq models. Figures |3] and H] show the time evolution in the HR 
diagram, and in the central density-temperature plane, respectively. In Figure HJ spikes are 
shown at Pc ~ 10^ g/cm^ for both 11.0 and 11.2 Mq models, representing off-center neon 
ignitions. In our calculation, the minimum initial mass for Ne ignition is 11.0 Mq and the 
CO core mass is 1.35 Mq. A star with a lar ger initial mass than this cri tical mass will form 
an Fe core and end up as a normal CCSN (INomoto fc Hashimoto! Il988l ). Since evolutional 
properties of both 10.4 and 10.6 Mq are similar to 10.8 Mq, here we mostly show results of 
a 10.8 Mq star, which has the largest core in these less massive stars and provides the most 
plausible progenitor model for an ECSN. 

The minimum ma ss for Ne ignition ~ 11 Mq is large compa red with other recent evo- 
lutionary calculations (jPoelarends et al.l 120081 : iPumo et al.l 120091 ). This is because our cal- 
culations do not take into account the effect of add itional mixing such as overshooting. If 
an exponentially decreasing diffusion (jHerwigl 120001 ) is treated, our calculation shows that 
this minimum mass is reduced by ~ 2 Mq with an overs hooting para meter of /over = 0.02. 
This behavior is quite consistent with the calculation by ISiesd (120071 ). In fact, inclusion of 
overshooting importantly changes the relation between initial masses and He core masses. 
While this severely affects an estimation of an initial mass range for ECSNe, effects on evo- 
lutionary results of progenitor calculation are much more less. Especially, due to a lack of 
observational constraints, most evolutionary calculati ons for massive stars only take into ac- 
count the overshooting before core C burning stages (IHirschi et al.ll2004j : iLimongi fc Chieffi 
20061 ) . In this case, our results on the progenitor evolution becomes fully consistent with 
these calculations. 
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3.1. Pre SAGB evolution 

In this section, we summarize the evolutional resuhs of a 10.8 star from its ZAMS 
phase to the completion of the dredge-out, or mergence of convective regions which surrounds 
the ONe core. Figure [5] shows the time evolution of convective regions. 

3.1.1. The hydrogen burning stage 

The duration of core hydrogen burning is 1.74x 10'' yr. In this stage, convection develops 
in the central region due to the luminosity generated by the CNO-cycle. Core hydrogen 
burning continuously shifts to shell hydrogen burning as the central hydrogen burns out, 
while the H-depleted core contracts and is heated up by the release of the gravo-thermal 
energy. The increasing luminosity in the outer shell burning region expands and cools the 
envelope. As the opacity increases, the convectively unstable region emerges at the surface. 
The base of the convective region extends inward, then the star becomes a red giant. When 
the first dredge-up episode occurs, ^^N, the second main product of the CNO-cycle, is dredged 
up to the surface. Surface composition of these isotopes at the end of this episode are 
summarized in Table [TJ 

3.1.2. The helium burning stage 

When the central density and temperature reach log pc = 3.50 and log Tc = 8.16, 
core helium burning takes place. The luminosity of the shell hydrogen burning decreases, 
and the base of the convective envelope retreats outward in mass. The whole region of 
the hydrogen envelope becomes convectively stable, and the star enters into a blue-loop 
in the HR diagram. The core helium burning phase continues 2.76x10^ yr, followed by 
shell helium burning after core helium depletion. Since the large luminosity by shell helium 
burning expands and cools the hydrogen burning layer, shell hydrogen burning dies out. This 
luminosity re-heats the envelope and convection arises in the envelope. The star becomes an 
AGB star, in which the partially-degenerate CO core has formed, surrounded by the helium 
burning shell. At the center of the core, the ratio of mass fractions of oxygen and carbon 
becomes A(C)/X(0) = 0.5728. 
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3.1.3. The carbon burning stage 

In the partially-degenerate CO core, an off-center C flashes take place. In our calcula- 
tion, the flrst two flashes (Ci & C2 in Fig. [6]) arise near the center and last seven flashes burn 
outward (from C3 to Cg in Fig. |6]). These C flashes transform core carbon into neon and 
other intermediate mass nuclei. After the end of the sixth C burning, the second dredge-up 
reduces the mass of the helium layer, and the extent of the second dredge- up is smaller tha n 



lower massive stars. These results are consistent with other calculations ( ISiessI l2006l . 120071 ). 
Figure [7| shows the evolution of the core structure in terms of /5e, the ratio of electron pres- 
sure to the total pressure, temperature, and density at flve different stages: disappearance of 
convective core helium burning (la), commencement of the flrst shell carbon burning (lb), 
ignition at the center of the core (Ic), ignition of the eighth C burning (Id), and occurrence 
of the dredge-out (le). 

In advance of the off-center C ignition, the gravo-thermal heating increases the tem- 
perature of the contracting CO core (Fig. [TJ la). When the central density reaches log 
Pc = 6.28 and the maximum temperature in the core reaches log Tmax = 8.81 at Mr = 0.05 
Mq, off-center carbon ignition takes place (Fig. [TJ lb). This is because relatively high tem- 
perature activates the cooling process via neutrino emission in the CO core. The neutrino 
cooling efficiently removes local heat and suppresses temperature increase. The efficiency 
becomes greater with higher density, so for the inner region of the core, the cooling becomes 
more effective. However, the inner region of the core is more degenerate and thus harder to 
contract. On the other hand, still mildly degenerated outer region liberates gravo-thermal 
energy and increases its temperature, owing to the contraction by neutrino cooling. As a 
result, an inverse temperature gradient takes place in a degenerate core. 

Since the flrst C burning expands and cools the central region, the burning front of 
this shell burning does not propagate inward and the central CO region remains. On the 
other hand, the second shell C burning, which takes place as the flrst burning dies out, 
propagates inward and the center of the core ignites at last (Fig. [TJ Ic). The surrounding 
region of the second C burning shell has lower mass fraction of ^^C owing to the flrst C 
burning. This weakens the second shell burning, resulting in a smaller expansion of the 



central region. This enables the second burning flame to propagate inward (ISiessll2006l ). The 
propagation is caused by heat conduction from the base of the nearly stationary flame, the 
mean propagation speed becomes 9.7 x 10^'^ cm/sec. Final abundances at the center of the 
core are shown in Table [2] in terms of mass fractions. 

During the following four C flashes (Cs-e), the core keeps contracting. At the ignition 
of the seventh C burning, the center of the core is supported only by perfectly degenerate 
electrons, and the central region clearly shows a temperature inversion (Fig. [TJ Id). On 
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the other hand, radiative pressure has a larger contribution to the total pressure both at 
the flame front (Mr=1.24 Mq) and at the edge of the core (Mr=1.35 Mq). This is due to 
the high temperature and low density at these regions. Figure [S] shows the evolution of the 
energy structure and convective region around the edge of the core. Core contraction heats 
the base of the surrounding helium burning shell (HeBS) and increases its luminosity. This 
He shell burning induces convection at the base of the HeBS (Fig. [HI 2a), and shell He 
burning keeps supporting the convection during the seventh and eighth C burning phases 
(Fig. El 2b & 2c). 

After the end of the eighth C burning, luminosity from the shell He burning takes its 
maximum value owing to the temperature rise at the base of the HeBS (Fig. |8l 2d). However, 
soon the luminosity decreases, and the convection in a helium layer starts to be supported 
by escaping luminosity of C burning in the core (Fig. |H1 2e), then the growing convective 
region in the helium layer merges with the outer convection in the hydro gen envelope . At 



the moment of convective mergence, or the dredge-out episode named by llberu (119971 ). the 
temperature and density show steep drops at the edge of the core from log T = 9.0 — )■ 6.8 and 
log p = 4.8 — )■ —2.4 in a narrow mass range of 5 x 10~^ Mq (Fig. [TJ le). Some hydrogen at 
the envelope are mixed into the base of the HeBS. The resulting hydrogen burning generates 
energy in a very high rate of ~ 10^^ erg/sec/g (Fig. [HI 2f). The ONe core mass, which 
is defined by the mass coordinate of the maximum energy generation of helium burning, is 
1.347 Mq at the end of the carbon burning stage. 

Since the dredge-out episode mixes CNO products with the hydrogen envelope as well 
as other dredge-up episodes, the surface composition changes after the convective mergence 
(Table [1]). After completion of the dredge-out, the ONe core is surrounded by thin He shell 
and diffuse H envelope. In order to calculate the further evolution of the ONe core, core 
mass growth by shell helium burning should be considered. Thermal pulses take place at 
the thin helium layer in the SAGB star as in the case for lower massive AGB stars, however, 
simulation for the thermal pulse phase requires high cost calcula t ion an d full calculation of 



the phase is difficult. Moreover, as shown by iPoelarends et al.l (120081 ) for the case of the 
dredge-out phenomenon, not only helium but also hydrogen may be mixed with the base of 
the envelope and get burned out with large amount of energy production. This requires a 
scheme which can solve both mixing and reaction simultaneously, however, our code is not 
equipped with such a scheme at present. So from the end of the dredge-out episode, we 
assume the constant shell helium burning and avoid the numerical difficulties on core growth 
(see § 2.5). The radiative- zero approximation is set to be the boundary condition to the 
edge of the core. Since the edge structure is extremely steep, this approximation would not 
affect the later core evolution, especially at the central region. 
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3.2. Evolution of a contracting ONe core 

In this section, we sliow ttiat tlie evolution of an contracting ONe core can be divided into 
four sub-phases in terms of driving mechanisms; neutrino cooling, core mass growth, electron 
capture on ^^Mg and ^*^Ne, and 0+Ne deflagration. The efliciency of these mechanisms are 
related to evolutionary timescales which are shown in Fig. [Hi The deflnitions of timescales 
are given by 



"^^-^ - rflnp,' ^™ - R.ore{L + Uy - rflnM,o,e' ^^^^ " rflnFe' " V GM,ore'^ ^ 

and they represents the timescale of core contraction, the Kelvin-Helmholtz timescale, the 
timescale of core mass growth, the timescale of electron capture, and the dynamical timescale, 
respectively. 



3.2.1. Contraction owing to neutrino cooling 

From the beginning of the contraction until the central density reaches log pc = 9.39, 
the contraction is caused by neutrino cooling. Figure |9] shows that the time scale of contrac- 
tion is close to the Kelvin-Helmholtz time scale, corresponding to the neutrino-cooling time 
scale. In this phase, the central temperature decreases as the entropy is radiated by ther- 
mally activated neutrino emission. As the central temperature decreases, the cooling rate 
of neutrino emission decreases as well. The timescale of density evolution becomes longer 
simultaneously. 



Although electron captures on ^''Al, ^^Mg and ^^Na proceeds in this stage (Fig. [T0|) . 
these reactions have just a minor effect on the density evolution. This is because, flrstly, 
mass fractions of these nuclei are so small that reduction of electron mole fraction does not 
induce contraction. Secondly, thermal contributions from these reactions are smaller than 
neutrino cooling. Since the energy production by th ese electron captures is small, convection 



and thus convective URCA cooling as described in iRitossa et al.l ( 119991 ) does not affect our 
calculation. 



3.2.2. Contraction owing to core mass growth 

As the core mass increases, the central pressure, which is required to support the core, 
and thus the central density increase as well. While the efficiency of neutrino cooling is 
suppressed by entropy reduction, the efficiency of core mass growth is independent from 
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core structure and becomes constant. After neutrino cooling becomes less effective than core 
mass growth, the timescale of core growth Xgrowth starts to limit the contraction timescale. 

The stationary core growth forces the central density to increase in a constant rate, 
and the production rate of gravo-thermal energy becomes constant. Owing to this constant 
heating, and also owing to less effective neutrino cooling, the central temperature increases 
proportionally to the central density. 



3.2.3. Contraction owing to electron captures 

After the central density exceeds log pc = 9.88, core contraction is driven by electron 
captures and thus, Teiec starts to rule the evolution time scale. The contraction time scale 
decreases with increasing density owing to the increasing electron capture rate. Soon core 
growth becomes negligible since Teiec becomes much smaller than Tgrowth- After that the core 
mass is frozen and is considered as the critical core mass for an ECSN, Mec- For a 10.8 Mq 
model, Mec is 1.366 Mq. 

Figure S] shows that the central temperature increases with the central density more 
steeply in this core growth stage. The increase of temperature is due to heating via electron 
captures and the main reaction sequence is ^^Mg — )■ ^^Na — )■ ^^Ne at that time. Electron 
capture on ^"^Mg forms an excited ^^Na* as a daughter nuclei. When this ^^Na* decays 
to the ground state, 7-ray accompanies and heats surroundings. Owing to this additional 
heating, the electron capture on ^^Mg becomes exothermic. Moreover, both ^^Na* and ^^Na 
at the ground state can capture another electron with lower threshold density than the one 
of ^"^Mg. Therefore, electron capture on ^"^Mg results in double-electron capture and releases 
large amount of heat. It is noteworthy that the second daughter nuclei ^^Ne becomes the 
dominant products of electron capture on ^^Mg (Fig. [T0|) . As a result, the amount of 
Vrad becomes sufficiently large compared to |V^, and the central region is fully mixed by 
convection. The growing convection supplies fresh ^^Mg to the center where electrons are 
quickly captured and the electron m ole fraction d e crease s in the whole convective region. 



In accordance with the discussion by iMiyaji et al.l ( 1l980l ). convective URCA by ^^Mg-^^Ne 



does not take place in this phase, though the reaction drives convection. This is because 
the threshold density of /3~-decay for the daughter nucleus ^"^Ne is ~ 10^ g/cm^ and is too 
small to occur in this phase. While minor electron capture nuclei, such as ^^Na and ^^Si, 
can induce convective URCA in this phase, the averaged energy loss by these reactions are 
reasonably smaller than heating by electron capture on ^"^Mg, and thus this process can be 
negligible. 
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After the central density reaches log pc = 10.3, ^''Ne starts to capture electrons (Fig. 
[TOl) . This electron capture keeps driving convection as well as ^^Mg. The time scale of the 
evolution becomes shorter than the one of convection, and partially mixing allows a gradient 
of X(^°Ne) to exist in the central region. However, since ^°Ne has a large mass fraction 
of ~0.4, until the commencement of Ne+0 deflagration, the fuel is not consumed and the 
reaction continues to heat the core. 



3.2.4. 0+Ne deflagration 

When the central temperature reaches log = 9.2, 0+Ne burning takes place at the 
center and the central temperature rises in a very short time scale. Since reaction rates 
of the nuclear burning highly depend on the temperature, thermal runaway takes place in 
the extremely degenerate region. While the temperature steeply increases at the moment 
of the ignition, the variance of the pressure and density becomes small. When the central 
temperature exceeds 5 x 10^ K, the NSE composition is assumed to be achieved. 

At the boundary of the NSE region and the surrounding ONe region, there must be a 
negative steep gradient of entropy as a result of the 0+Ne burning, and convection must 
exist. In our calculation, heat transportation by convection is solved (see § 2.6), and the 
resulting heating at the base of the ONe region becomes more efficient than the one of 
electron captures on neon and magnesium. The heating increases the temperature, soon the 
0+Ne burning ignites and the NSE is achieved. Top panels of Fig. [11] show the evolution 
of temperature distribution. The location of the steep temperature gradient is identical to 
the burning front and propagates outward. The propagation velocity, ~ 1.0 x 10^ km/sec, 
becomes smaller than the sound velocity. In other words, the NSE region extends outward 
in mass by the 0+Ne deflagration. 

In the NSE region, electron ca pture reactions on both free-protons and heavy nuclei 



take place (iJuodagalvis et al.ll2010l ). The reduction of the electron mole fraction affects 
the dynamical evolution of the core. Contrary to electron captures on ^^Mg and ^*^Ne, the 
electron capture in the NSE region becomes endothermic reaction. This is due to the global 
compositional change in the NSE composition; unstable neutron rich isotopes are preferable 
to exist in the lower environment, and thus reduction of Yg causes the free energy to be 
restored in terms of nuclear binding energy. Because of the cooling, the positive entropy 
gradient appears in the NSE region. The center of the core becomes convectively stable. 

The ONe core has been strongly bound by gravity and has a binding energy of —6.534 x 
gj,g g^j^^ total energy of —5.791 x 10^° erg at the commencement of 0+Ne ignition. 
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While energy injection by the 0+Ne burning is too small to disrupt the whole star, fast 
reduction of electrons in the NSE region accelerates core contraction. As a result of both 
energy injection and electron reduction, the contraction timescale of the core slowly decrease 
during propagation of the NSE region. After 2.36xl0~^sec from the ignition at the center of 
the core, the central density reaches log pc = 11.0, the deflagration front reaches Mnse = 0.12 
Mq and tnse = 1-51 x 10~^ R©, the binding energy and total energy of the core becomes 
-7.739 X 10^^ erg and -9.158 x 10^° erg, respectively. 



3.3. Parameter dependences on Mec 

In order to investigate the uniqueness of the critical core mass for ECSNe, Mec, in 
various situations, additional test calculations on ONe core contraction are done with dif- 
ferent settings. As an initial condition, ONe cores of 1.346, 1.332, and 1.288 Mq are taken 
to be formed in SAGB stellar models of 10.8, 10.6, and 10.4 Mq. For each ONe core, later 
contraction phase is calculated with three different rates of core mass growth of 1.0 x 10~^, 
1.0 X 10~^, and 1.0 x 10^'' Mq/jt and with or without the Coulomb correction for electron 
captures. Results are shown in Fig. [121 

Core masses at the end of calculations were taken as Mec for each model, however, 
calculations were stopped for several models owing to numerical difficulties during its electron 
capture phase. This is due to convective URCA processes of minor nuclei such as ^^Na and 
^^Si (see § 4.2 for discussion). For these models, we took core masses at the emergence 
of convection as Mec- Because the contraction timescale steeply decreases after electron 
capture on ^^Mg starts, and the central convection is driven by this reaction, this procedure 
would cause small errors on our results. For instance, the core growth after the emergence 
of central convection for Mini =10.8 Mq with standard settings is 0.002 Mq. 

The mean value of our results with the Coulomb correction is 1.367 Mq. We estimate 
the uncertainty of Mec of about ± 0.005, taking into account the error from the procedure 
of determination of Mec and the variances of the results. When the uncertainties which 
comes from the C burning phase are considered, the uncertainty should increase in total. 
We assume that the amount of this uncertainty is as large as one comes from later stages of 
evolution. Thus the uncertainty of our result is about ± 0.01 Mq in total. 
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4. Discussions and Conclusions 
4.1. The critical core mass for ECSNe 

In our calculation, M-^c is 1.367 Mc^ wit h a relatively small error of ±0.01 Mq, and is 



consistent with the result by iNomotd (Il987f ). M^c = 1-375 Mq. This result supports the 
idea that the critical core mass for ECSNe exists uniquely. I.e., in a critical object with the 
Chandrasekhar mass with given compositions, the central density can be represented only 
by the core mass, pc = Pc(^core)- Therefore, the critical mass for ECSNe should be uniquely 
determined by solving the relation of /324Mg = Pc{M^c), where pz^ug is the threshold density 
for electron capture on ^^Mg. 

Our result shows that inclusion of the Coulomb correction increases pziug by ~10% and 
increases Mec as well. Also both faster and slower core growth may increase Mec- How- 
ever, these effects are smaller than the numerical uncertainties and the detailed quantitative 
analysis is difficult. This is because, with these critical objects, even a tiny difference in 
a physical quantity can cause an important change in the results. For instance, change in 
temperature distribution will cause a change in Mec- K is difficult to uniquely determine 
the distribution of temperature, because the temperature of degenerate electrons is almost 
independent of the dynamical structure. Thus the temperature should be determined only 
by nuclear reactions, which relatively have larger uncertainties. The critical masses for 10.6 
Mq models show larger values than the other two initial masses. This may be due to the 
different distribution of temperature and chemical compositions. 



4.2. Convective URCA by minor electron capture nuclei 

When electron capture on ^"^Mg drives convection, convective URCA process is induced 
by minor electron capture nuclei such as ^^Na and ^^Si. In a central convective region, 
daughter nuclei, which are formed at the center of the core, are mixed with outer less dense 
region. Sometimes the density at the outer edge of the convective region is much less than 
the threshold densities for /3~-decays of these nuclei. In such cases, reaction rates of 
decays become significantly large, and this results large energy loss only at the edge of the 
convection. Some of our models are stopped by numerical difficulties caused by this energy 
loss. 

However, this behavior seems unphysical, and more slow reactions will occur in reality. 
This is because such a large energy loss should immediately stop the convective mixing 
around the region, and mass fractions of these nuclei are so small that durations of these 
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reactions should be limited in a non-convective environment. Therefore, if small fractions 
of unstable /3~-decay nuclei are mixed into less dense regions, the resulting endothermic 
reaction will stop immediately and have minor effects on the whole evolution. We expect 
that the convective URCA process during the electron capture phase will not affect the later 
evolution of an ONe core as shown in the 10.8 M© model. To proof this statement, it is 
required to solve the coupling of convective mixing, nuclear reactions, and a stellar structure 
simultaneously. 



4.3. Deflagration and core-collapse 



Figure |9] shows that the contraction timescale Tcon is longer than the dynamical timescale 
T"dyn even at the end of calculation. Such quasi-static contraction will be altered into dynam- 
ical collapse during the deflagration phase, and thus, the core-collapse in an ECSN is caused 
by electron captures in the NSE region after the initiation of 0+Ne deflagration, rather than 
by electron captures on ^^Mg and ^°Ne. In order to provide a plausible progenitor model for 
an ECSN, the deflagration phase should be considered carefully because propagation of the 
deflagration front and both energy generation and electron captures in the NSE region can 
affect the core evolution importantly. 

The obtained deflagration velocity 1.6x10^ km/sec at the e nd of the calc u latiori (pc = 
10^^ g/cm^) is consistent with a velocity ~ 10'^ km/sec as in iMiyaji et al.l ( 1l980l ). On 
the other hand, th e extension of the NSE reg i on, 0. 12 Mq, is smaller than 0.354 M0 by 
Miyaji et al.l (1l980[ l and ~0.3 M© by iNomotd (1l987l ). in which the same assumptions for 
convective energy transportation and electron capture rates are taken. This i s due to a 
different treatm ent of electron captures in the NSE region. Both iMiyaji et al.l (Il980l ) and 
Nomotd ( 119871 ) assumed that only free-protons can capture electrons in the NSE compo- 
sitions, and thus electron capture rates at low Ye environment with extremely low proton 
fractions become much smaller than ours , which takes into account electron captures by neu- 
tron rich nuclei ( jJuodagalvis et al.ll2010( ). In order to confirm the effect of electron captures 
on neutron-rich nuclei, we calculate the deflagration phase taking into account electron cap- 
ture only by free-protons in NSE compositions. The resulting propagation of the NSE front 
is shown in Fig. [131 Comparing with the case of full electron captures, contraction becomes 
slower owing to smaller reduction of electrons. Then the duration of contraction becomes 
longer, and the flame front propagates far from the center at the end of the calculation even 
though the flame front has a smaller mean velocity of 8.4 x 10^ k m/ sec. Since the explosion as 
an ECSN is triggered by neutrino heating f lKitaura et al.ll2006l ) and the efficiency is affected 
by the density profile, a central concentrated density profile in our calculation would affect 
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its explosion and the inclusion of electron captures by neutron-rich nuclei will be important. 



4.4. Constraints by mass loss 

It is considered that not all of SAGB stars end up as ECSNe because of intense mass loss 
during SAGB phase. As the core increases its mass, the outer envelope is lost by this mass 
loss. This limits the duration of SAGB phase, and if the core mass is less than Mec at the 
end of SAGB phase, the star ends up as an ONe WD. Therefore the ratio C, = Mcore/|M;nv| 
is important to determine stellar fates. 

The minimum ^, with which a star becomes an ECSN just before the star completely 
loses its envelope, can be defined as 

_ Mec - MZe 

where M™J.g and M^l represent an ONe core mass and a total mass at the beginning of SAGB 
phase. When an actual ^ for a SAGB star is larger than ^crit, the star ends up as an ECSN, 
and vice versa. If core growth rates are specified by certain simulations, this prescription 
becomes identical to the definition of the critical mass loss rate |Menv,crit|, thus 

|Menv,crit| = ^cri\"^core; (28) 

and the star with |Menv| < |Menv,crit| can become an ECSN. 

^crit becomes 2.41x10-3, 4.01x10^3, and 9.25 xlO'^ for 10.8, 10.6, and 10.4 Mq models 
respectively. Because of large uncertainties on both Mmrp. and Mt^t, current estimate of ^ 



is highly uncertain. Here, according to the results by iPoelarends et al.l (120081 ). we limit the 
range of ^ as 2.54x10-^-1.80x10-^. This estimate includes wide ranges of mass loss rates, 
initial masses, and thus mass growth rates, and the minimum value will be too small to be 
applied to the most massive SAGB stars in our calculation. Even from this crude estimate, 
our model of a 10.8 Mq star is plausible for the progenitor of an ECSN. This conclusion will 
be robust if the inclusion of overshooting reduces the total mass by ~2 Mq and rises ^crit to 
3.13x10-3. 



4.5. Conclusions 

Stellar evolution for the most massive SAGB stars with solar abundances were calculated 
in this work, and we provided progenitor models for ECSNe with detailed evolutionary 
calculation for the first time. 
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In such a star, nuclear burnings of hydrogen, hehum, and carbon take place step by 
step, then a core mainly made of oxygen and neon forms. The ONe core, which is supported 
by degenerate electron gas, contracts due to neutrino cooling, core mass growth by the sur- 
rounding shell helium burning, and reduction of electron mole fraction by electron capture 
reactions. When the central temperature rises high enough to ignite oxygen, 0+Ne defla- 
gration takes place and NSE is achieved in the region. Although the nuclear burning heats 
the region, the released energy is too small to explode the highly gravitationally bound core. 
As electron capture reactions on neutron-rich nuclei as well as on free protons accelerate 
the core contraction, the deflagration front propagates outward. The core will continue to 
contract up to the formation of a proto-neutron star. The fate will be a weak type II SN. 

In our results, the critical mass for ECSNe, M-^c, is 1.367 Mq with a relatively small 
uncertainty of about ±0.01 Mq, and we showed the uniqueness of the value under various 
initial core masses and various core growth rates. The uncertainty may come from the 
computational difficulties on the physically critical object as discussed in § 4.1. Inclusion of 
the Coulomb correction may increases Mec, however, the effect is smaller than the numerical 
uncertainties. 

Inclusion of minor intermediate mass nuclei such as ^^Si, ^^Na, ^^Mg, and ^^Al does not 
affect the core evolution and thus Mec- Though convective URCA process by these nuclei 
sometimes cause numerical difficulties in an electron capture phase, this process will have 
minor effect on the whole evolution. 

Since the ONe core keeps contracting quasi- statically even during the deflagration phase, 
the core evolution is importantly affected by propagation of the deflagration front and elec- 
tron capture in the NSE region. We showed that the assumption of electron capture only by 
free-protons leads slower contraction and electron captures on neutron-rich nuclei should be 
treated. 

Owing to intense mass loss during SAGB phase, not all of our models may end up as 
ECSNe. Accurate estimates for the core growth rate and the mass loss rate are difficult and 
further investigations are still needed on t his topic. Howeve r , even under the most severe 



condition which is taken from the work by iPoelarends et al.l (|2008[ ). our 10.8 Mq model is 
plausible for a progenitor model of an ECSN. 
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Table 1: Surface abundances at three different stages of ZAMS, the 1st dredge-up, and the 
dredge-out. 
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Fig. 1. — Nuclei with proton number less than 14 that are included in our calculation. 
Eighteen blue nuclei are newly included and 300 in total are taken into consideration, while 
red nuclei had been taken into account in our previous paper (Y&Ull). 
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Fig. 2. — The NSE composition on (N, Z) plane in the cases of (temperature [10^ K], density 
[g/cm^], electron mole fraction) = (10.0, 10^, 0.49) and (10.0, 10^^, 0.27). Blue region shows 
the region of considered 3091 nuclei and red boxes represent stable nuclei. Color contours 
show the mass fraction distribution. 
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Fig. 3. — Evolution of ~11 Mq models in the HR diagram. 
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Fig. 4. — Evolution of ~11 Mq models in the central density-temperature plane. 
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Time until the dredge — out (year) 

Fig. 5. — The time evolution of convective regions of a 10.8 Mq model until completion of 
the dredge-out. Convective regions are shown in green hatched region in the figure, and log 
scaled net nuclear energy generation is shown in color. 




Time until the dredge — out (year) 
Fig. 6. — Same as Fig. |5]but expanded on a off-center C burning phase. 
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Fig. 7. — Structure evolution in terms of electron pressure fraction (top), temperature (mid- 
dle), and density (bottom) as functions of mass coordinate at five different stages: disap- 
pearance of convective core helium burning (la), commencement of the first shell carbon 
burning (lb), ignition at the center of the core (Ic), ignition of the eighth C burning (Id), 
and occurrence of the dredge-out (le). Points show the locations of the base of the HeBS for 
(Id) and (le). 



Table 2: Central abundance at the end of core carbon burning phase. 
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Fig. 8. — Profiles of tfie energy generation rates and luminosity witli mass coordinate at 
different stages: occurrence of convection on HeBS (2a), ignition of tlie seventli C burning 
(2b), ignition of tlie eigfitfi C burning (2c), shell He burning at its maximum (2d), 6.84 x 
10 yr after (2d) (2e), and occurrence of the dredge-out (2f). Different line types are the 
nuclear energy generation (red, solid line), the neutrino energy loss (green, dashed line) the 
gravo-thermal energy release (blue, dotted line), and the luminosity (cyan, dash-dotted line) 
respectively. Gray regions represent convective regions. 
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Fig. 9. — Evolution of the timescales with density at the stellar center. Different lines 
correspond to different timescales: timescale of core contraction Tcon (red, solid line), Kelvin- 
Helmholtz timescale tkh (green, long dashed line), timescale of core mass growth Xgrowth 
(blue, short dashed line), timescale of electron capture Xeiec (magenta, dotted line), and 
dynamical timescale Xdyn (cyan, dash-dotted line). Background colors shows different sub- 
phases: the neutrino cooling phase (log Pc<9.39), the core growth phase (9.39<log Pc<9.88), 
the electron capture phase (9.88<log pc< 10.39), and the deflagration phase (10.39<log pc). 
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Fig. 10. — Evolution of mass fractions of electron capture nuclei at the stellar center with 
central density. For lower density, mass fractions of ^^Mg, ^^Na (red & green, solid line), 
^^Na, ^^Ne (red & green, dashed line), and ^^Mg (orange, dotted line) are shown. For higher 
density, mass fractions of ^^Mg, ^^Na, ^"^Ne (orange, green, & blue, dotted line) and ^°Ne, 
20p^ 20q (^orange, green, & blue, dash-dotted line) are shown. 
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Fig. 11. — Profiles of temperature, density and electron mole fraction during 0+Ne defla- 
gration are shown as a function of both mass and radius coordinates. These profiles are 
taken at 1.13 x lO'^, 5.40 x IQ-^, 9.92 x 10''^, 1.49 x 10-\ 2.00 x 10-\ and 2.34 x 10"^ sec 
after the ignition at the center of the core. 




Fig. 12. — The critical core masses for ECSNe calculated with different settings. Different 
colors show different core growth rates and different point types show different treatments 
of the Coulomb correction. 
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Fig. 13. — Same as Fig. [11] with electron capture rate only by free-protons. These profiles 
are taken at 1.15 x 10"^ 5.56 x 10"^, 1.00 x lO'^, 2.00 x 10-\ 3.00 x 10-\ and 3.58 x 10"^ 
sec after the O+Ne ignition at the center of the core. 



